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Abstract 

In this article we have calculated the structure properties of a strange quark star in static model 
in the presence of a strong magnetic field using MIT bag model with a density dependent bag 
constant. To parameterize the density dependence of bag constant, we have used our results 
for the lowest order constrained variational calculation of the asymmetric nuclear matter. By 
calculating the equation of state of strange quark matter, we have shown that the pressure of this 
system increases by increasing both density and magnetic field. Finally, we have investigated the 
effect of density dependence of bag constant on the structure properties of strange quark star. 



* Corresponding author. E-mail: bordbar@physics.susc.ac.ir 



1 



I. INTRODUCTION 



The core of a neutron star has been formed from nuclear matter, composed of the neu- 
trons, protons, electrons (for assurance negati on of electric charge) and other particles like 
pions, mesons and etc (Lattimer & Prakash 2004J). It is known that nuclear matter is 
meta-stable, that after releasing a lot of energy converts into strange quark matter (SQM) 
to achieve stability. This quark matter is the most stable state of matter that has been 
known until now. Thus, there is a new class of compact stars that come from the collap se 
of neutron stars, and are more stable compared to the neutron stars (Farhi & Jaffe 1984] ) . 
The best candidates for this conve rsion are t he ne utron st ars wi th masses of 1.5 — 1.8 M Q 
and quick spin (Drake et al. j2002j |: Li et al. [19991 ] : Weber [20051]). 

The collapse of a neutron star may lead to a strange quark star (SQS) or a hybrid star. 
Also under special conditions, an SQS may be directly born from the core collapse of a type 
II supernova. An SQS, from its center to surface is ma de from SQM, and on its surface 
may exist a nuclear layer (Glendenning & We ber 19921 ]) . Hybrid stars are the ones with 



cores composed of SQM (Bhattacharyya et al. [2006| ) . Here, we just consider the structure 
properties of SQS. 

The mass and density of an SQS is between the mass and density of a neutron star and 
that of a black hole. The mass-radius relation for an SQS is as M oc B? which is different 
from that of a neutron star. This star does not have the minimum mass. For an SQS 
with 1M < M < 2M Q , the radius is about 10 km (Farhi & Olinto [19861 ] : Shapiro & 
Tenkolsky [19831 ] V 

Rec ent ob servations indicate that the object SWIFT J1749. 4-2807 may be an SQS (Yu 
& Xu 201CH ] ) . The given results by Chandra observations also show that the objects RX 
J185635-3754 and 3C58 may be bare strange stars (Prakash et al. 20031 ] ). It is known 
that the compact objects such as neutron star, pulsars, magnetars and strange quark stars, 
are under the influence of strong magnetic fields whic h are typically abou t 10 15 — 10 19 G 



(Kouveliotouet al. [1998) . [1999) : Haensel et al. [2007] : Glendenning [2000| : Weber [20071 ]: 



Camenzind 



20071 ] ) . Therefore, in astrophysics, it is of special interest to study the effect of a 



strong magnetic field on the properties of SQM. We note that in the presence of a magnetic 
field, the conversion of neutron stars to bare quark stars can not ta ke place unless the value 
of the magnetic field exceeds 10 20 G (Ghosh & Chakrabarty 20011 ] ). 
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In recent years, we have calculated the maximum gravitational mass a nd oth er structure 
properties of a neutron star with a q uark c ore at zero (Bordbar et al. 2006] ) and finite 
temperatures (Yazdizadeh & Bordbar [201 We have also computed the structure prop- 
erties of SQ S at zero temperature (Bordbar et al. 20091 ] ) and finite temperature (Bordbar 



et al. |201l| ). We have also calculated the structure of a mag netize d SQS using MIT bag 
model with a fixed bag constant (90 y^jr) (Bordbar & Peyvand 201 In the present work, 
we investigate the effect of density dependence of bag constant on the structure of an SQS 



in the presence of strong magnetic field. 



II. COMPUTATION OF STRANGE QUARK MATTER EQUATION OF STATE 
IN THE PRESENCE OF MAGNETIC FIELD 



The equation of state (EOS) of strange quark matter (SQM) plays an important role 
for determining the structure of stars at high densities. To obtain EOS of SQM, there are 
different models based on Quantum Chromodynamics (QCD). At present, it isn't possible to 
achieve an exact EOS of SQM by primary principles of QCD. Therefore, scientists have tried 
to find approximate meth ods by combini ng; th e basic features of QC D, for exampl e, M 



bag model (Chodos et al. [1974 : Weber [2003]; Peshier et al. [2000j; Alfordet al. [2005) ). 



N JL model (Rehberg et al. [l996| ; Hanauske et al. 



200 



et al. 
Jaffe 



1984|) 



20061 ). and perturbative QCD model (Baluni 



; Ruster & Rischke 



1978) : Fraga et al. 



2004 1; Menezes 



20011 ]: Farhi & 



In MIT bag model, the quarks in the bag are considered as a free Fermi gas, and the energy 
per volume for SQM is equal to the kinetic energy of the free quarks plus a bag constant (B) 
(Chodos et al. 1974J). The bag constant B can be interpreted as the difference between the 
energy densities of the noninteracting quarks and the interacting ones. Dynamically, its role 
is as the pressure that keeps the quark gas in constant density and potential. In the initial 
MIT bag model, different values such as 55 and 90 j^f are considered for the bag constant. 
As we know, the density of SQM increases from surface to the core of SQS, therefore using 
a density dependent bag constant instead of a fixed bag constant is more suitable. 
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A. Density dependent bag constant 



The analysis of the experimental data achieved at CERN shows that the quark-hadron 
transition happens at a de nsity about seven times t he no rmal nuclear mat ter en ergy density 



;i56 MeVfm- 3 ) (Heinz [200 1) ; Heinz & Jacobs 



2000| ; Farhi & Jaffe 



1984J ). However 



theoretically, for no density-independe nt val ue of bag constant the hadron to quark matter 
transition takes place (Burgio et al. 20021 ] ). Therefore, it is essential to use a density 
dependent bag co nstant . Recently, a den sity d ependent form h as bee n also considere d for B 
(Adami & Brown [19931 ]: Jin & Jenning [19971 ]: Blaschke et al. |l999j ]: Burgio et al. |2002j). 
The density dependence of B is highly model dependent. According to the hypothesis of a 
constant energy density along the transition line, Burgio et al. tried to determine a range of 
possi ble val ues for B by exploiting the experimental data obtained at the CERN SPS (Burgio 



et al. 



2002]). By assumption that the transition to quark-gluon plasma is determined by the 



value of the energy density only, they estimated the value of bag constant and its possible 
density dependence. They attempted to provide effective parameterizations for this density 
dependence, trying to cover a wide range by considering some extreme choices in such a 
way that at asymptotic densities, the bag constant has some finite value. They employed a 
Gaussian form as follows 



B(p) = B oo + (B -B oo )e 



Mp/po) 2 



(1) 



400 In the above 



The parameter £> = B{p = 0) is constant and equal to B — 
equation, 7 is a numerical parameter which is usually equal to po ~ 0.17 fm~ 3 , the normal 
nuclear matter density. B^ depends only on the free parameter Bq. 

The value of the bag constant (£>) should be compatible with experimental data . Th e 
experimental res ults a t CERN-SPS c onfirm s a proton fraction x p t = 0.4 (Heinz 20011 ]; 
Heinz & Jacobs jsOOO | : Burgio et al. 2002 ] ). Therefore, we use the equation of state of 
asymmetric nuclear matter to evaluate B^. We use the lowest order constrained variational 
(LOCV) many-body method based on the cluster expansion of the energy for calc ulatin g 
;he eq uation of state of asym metric nuclear matter as fol lows (B o dbar fc Modarr es 1997], 



1998) : Moda rres fc Bordbar )l998| : Bordbar & Bigdeli |2007aj, |2007bJ, |2008a), |2008bj; 
Bigdeli et al. [2003 ]; Bigdeli et al. |201o| ). 

The asymmetric nuclear matter is defined as a system consisting of Z protons (pt) and 
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N neutrons (nt) with the total number density p = p pt + p n t and proton fraction x pt = ^y, 
where p pt and p nt are the number densities of protons and neutrons, respectively. For this 
system, we consider a trial wave function of the form, 

i> = Ft, (2) 

where is the slater determinant of the single-particle wave functions, and F is the A-body 
correlation operator (A = Z + N) which is given by 

F = SUf(ij). (3) 

i>j 

In the above equation, S is a symmetrizing operator. 

For the asymmetric nuclear matter, the energy per nucleon up to the two-body term in 
the cluster expansion is as follows 

The one-body energy, E\, is 

i=l ki 

where labels 1 and 2 are used for proton and neutron respectively, and ki is the momentum 
of particle i. The two-body energy, E2, is given by 

£2=^£<u|V(12)|ij-jz>, (6) 

where 

V(12) = "^[/(12), [V? 2 , /(12)]] + /(12)V(12)/(12). (7) 

In Eq. (j7j), /(12) and V(12) are the two-body correlation and nucleon-nucleon potential, 
respectively. In o ur calc ul ations, we use UVu + TNI nucleon-nucleon potential (Lagaris 



& Pandharipande [1981a| . [1981bl |). We minimize the two-body energy with respect to the 
variations in the correlation functions subject to the normalization constraint. From the 
minimization of the two-body energy, we get a set of differential equations. By numerically 
solving these differential equations, we can calculate the correlation functions. The two-body 
energy is obtained using these correlation functions and then we can calculate the energy of 
asymmetric nuclear matter. The proc edure of these calculations has been fully discussed in 



reference (Bordbar & Modarres 



1998^ . 



The experime ntal r esults at CERN -SPS confirm a proton fraction x p t = 0.4 (Heinz 20011 ] : 



Heinz & Jacobs 



20001 ]; Burgio et al. [20021 ] ). Therefore, to calculate B^, we use our results 



of the above formalism for the asymmetric nuclear matter characterized by a proton fraction 
Xpt = 0.4. By assuming that the hadron-quark transitio n take s place at the energy density 
equal to 1100 MeVfm' 3 (Heinz [20011 ]; Burgio et al. [20021 ] ). we find that the baryonic 
density of nuclear matter corresponding to this value of the energy density is ps = 0.98 /m~ 3 
(transition density). At densities lower than this value, the energy density of SQM is higher 
than that of the nuclear matter. With increasing the baryonic density, these two energy 
densities become equal at the transition density, and above this value, the nuclear matter 
energy density remains always higher. Later, we determine B^ = 8.99 by putting the 
energy density of SQM and that of the nuclear matter equal to each other. 



B. Energy density calculation of strange quark matter in the presence of magnetic 
field 



We consider SQM composed of u, d and s quarks with spin up (+) and down (-). We 
note the number density of quark i wit 
introduce the polarization parameter £j by 



denote the number density of quark i with spin up by p^ and spin down by p\ ' . We 
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pi 



(8) 



where < & < 1 and pi = p^ + p\ . Under the conditions of beta-equilibrium and charge 
neutrality, we get the following relation for the number density, 

P = Pu = Pd = Ps, (9) 

where p is the total baryonic density of the system. 

Within the MIT bag model, the total energy of SQM in the presence of magnetic field 
(B) can be written as 

E tot = E K + B + E M , (10) 

where Em is the contribution of magnetic energy, B is the bag constant (in this article, we 
use a density dependent bag constant (Eq. ([T])), and E K is the total kinetic energy of SQM. 
The total kinetic energy of SQM is as follows, 



E 



K — £^ 

i=u,d,. 



(11) 



where E\ is the kinetic energy of quark i, 

= E E \/ m i c4 + h 2 k(p) 2 c 2 (12) 
p=± k(p) 

we ignore the masses of quark u and d, while we assume m s = 150 MeV for quark s. After 
performing some algebra, supposing that £ s = £ u = = £, we obtain the following relation 
for the total kinetic energy density (e# = -^) of SQM, 



^ = -^-3 J2 [44 P) 4 P) (2^ 2 4 P)2 c 2 + my) - mfrftn( ^ + ^ /c )] 

107T Ai „ i_ C TTloC 



where 



and 



+ ^^p 4/3 [(l + 4/3 + (l-0 4/3 ] (13) 



fcW = (tt 2 p) 1/3 (1 ± 1/3 , (14) 



4 ± ) = (^4 ± ) 2 c 2 + m 2 c 4 ) i /2> (15) 



For SQM, the contribution of magnetic energy is as Em = —M.B. If we assume the magnetic 
field is along the z direction, the contribution of the magnetic energy of SQM is given by 

e m = - E m ? )b > ( 16 ) 

i=u,d,s 

where Mf> is the magnetization of the system corresponding to particle i which is given by 

M« = Niinti. (17) 



In the above equation, iVj and fa are the number and magnetic moment of particle i, respec- 

V 



tively. By some simplification, the contribution of the magnetic energy density (em = 4/ 1 ) 



of SQM can be obtained as follows 

£ M = ~ E PWiiiB- (18) 

i=u,d,s 

Using the above equation and p = p u = Pd = Ps an d with the assumption that £ = £ n = 
£d = £ s , we have 

£m = -(pB^p s + pB£p u + pB^p d ). (19) 

Now, w e take advantage of numerical values of the magnetic moment for quarks ( Wong 



Samuel 



19901 ) : 



p s = —0.581 Pn,Pu = 1-852 pN,Pd = —0.972 p^- 
Using Eq. (JT9~|) and above values, we conclude that 

e M = -0.299p£p N B, (20) 

where pn = 5.05 x 10~ 27 (J/T) is nuclear magneton. 

C. The results of energy of strange quark matter in the presence of magnetic field 

We have calculated the properties of strange quark matter in the presence of magnetic 
field with the density dependent bag constant (Eq. ([I])). Our results are as follows. 

Our results for the total energy density {e to t — ^r 1 ) of strange quark matter (SQM) in 
the presence of magnetic field have been plotted versus the polarization parameter in Fig. 
CD for various densities at B = 5 x 10 18 G. We see that there is a minimum point in the 
energy curve for each density which shows a meta-stable state for this system. As the density 
increases, the minimum point of energy shifts to the lower values of the polarization, and 
finally it disappears at high densities in which the system becomes nearly unpolarized. 

In Fig. [21 we have plotted the polarization parameter corresponding to the minimum 
point of energy versus density for two magnetic fields B = 5 x 10 18 G and B = 5 x 10 19 G. 
We can see that the value of £ decreases by increasing the density, and it becomes nearly 
zero at high densities. We have also drown the polarization parameter as a function of 
the magnetic field at different densities in Fig. [3j As this figure shows, the polarization 
parameter increases by increasing the magnetic field for all densities. 

For the strange quark matter, our results for the total energy density at B = 5 x 10 18 G, 
which calculated with the density dependent bag constant, have been shown as a func tion o f 
density in Fig. SJ The results for B = 90 W at B = 5 x 10 18 G (Bordbar & Peyvand [201 j ) 
are also given for comparison. It can be observed that the total energy density has an 
increasing rate by increasing the density. Also, it can be found that for p greater (lower) 
than about 0.6 fm~ 3 , the energy of SQM with the density dependent bag constant is lower 
(greater) than that with the fixed bag constant. From Fig. HJ it is seen that for p < 0.6 fm~ 3 , 
the increasing of energy has a slow slope, whereas for p > 0.6 fm~ 3 this increasing is 
accomplished with a more quick slope. 

Fig. [5] shows the phase diagram for the strange quark matter. We can see that as the 
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density increases, the magnetic field grows monotonically. It explicitly means that at higher 
densities, the ferromagnetic phase transition occurs at higher values of the magnetic field. 



D. The equation of state for strange quark matter in the presence of magnetic 
field 

In this section, we calculate the equation of state (EoS) of strange quark matter (SQM) 
in the presence of magnetic field with density dependent bag constant. Generally, we can 
calculate EoS using the following relation, 

P{p) = P^f - etot, (21) 

where P is the pressure and e tot is the energy density, which in the presence of magnetic 
field, is obtained from Eq. ( [TO]) . In Fig. El we have compared our results for EoS of 
SQM at different magnetic fields. This shows that for all magnetic fields, by increasing the 
density, pressure has an increasing rate. Also, we can see that with increasing magnetic 
field, the pressure increases. In Fig. we have drawn EOS of SQM by density dependent 
bag constant at B = 5 x 10 18 G. The results for B = 90 Ms£ at B = 5 x 10 18 G (Bordbar 



/m 3 

& Peyvand 20111 ]) are also given for comparison. This figure indicates that for p greater 
than about 0.52 fm~ 3 , when the bag constant is density dependent, the pressure of SQM is 
greater than that of the density independent case. 

III. STRUCTURE OF STRANGE QUARK STAR 

Quark stars are relativistic objects, therefore we used the general relativity for calculation 
of their structures. Since most of the massive general relativistic objects have some forms 
of rotation (very rapid in the case of pulsars). In this calculations, we are interested in the 
investigation of the strong magnetic field effects on the structure of a static strange quark 
star. Using the the equation of state (EoS) of strange quark matter (SQM), we can obtain 
the structure of these stars by numerically integrating the general relativistic equations 
of hydrost atic e quilibrium, the Tolman-Oppenheimer-Volkoff (TOV) equation (Shapiro & 



Teukolski [1983j ) 

,n ^/ £ ( r W r ) Wl l j , 

c 2 e(r)''^~ L ~~ r m(r)c 



dP G(^^)(l + 5fl)(l + ^M' 



dr (i _ 2G Mr) ^ 
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—— = 4vrr e(r). 
dr 



(22) 



In the above equations, P is pressure and e(r) is the energy density, G is the gravitational 
constant and m(r) is the mass inside radius r which is calculated as follows 



mir 



4nr' 2 e(r')dr'. 



(23) 



Now, by selecting a central energy density e c , under the boundary conditions P(0) = P c and 
m(0) = 0, we integra te TO V equations outwards to a radius r = R at which P vanishes 
(Shapiro & Teukolski |1983| ). 

In this section, we calculate the structure of the SQS with density dependent bag constant 
in the presence of a magnetic field. We should note that a strong magnetic field changes 
the spherical symmetry of the sy stem a nd for magnetic fields les s tha n 10 19 G, this effect 
is ignorable (Felipe & Martinez 20091 ] ; Gonzalez Felipe et al. 20111 ]) . Considering the 
anisotropy of the strange quark matter pressures in the presence of magnetic field, it has 



been shown that for vanishing AMM (anomalous magnetic moments), the perpendicu 



ar 
])• 



component of the pressure P± goes to zero at about 2 x 10 19 G (Gonzalez Felipe et al. [20081 
Thus in the case of SQM, for B < 10 19 G the anisotropy in the pressures is relatively small, 
i.e, P± = P\\. 

In Fig. [HI we have drawn the gravitational mass versus the central density (e c ) for an SQS 
in the magnetic fields B = and 5 x 10 18 G. We see that as the central density increases, 
the gravitational mass of an SQS increases, and finally it reaches a limiting value which is 
called maximum gravitational mass. Fig. [8] shows that by presenting the magnetic field, the 
gravitational mass decreases. The results for B = 90 *ff$ at B = 5 x 10 18 G (Bordbar & 



Peyvand 



20111 ] ) are also given in Fig. [8] for the sake of comparison. This indicates that the 



density dependence of bag constant leads to substantially higher values for the gravitational 
mass of SQS. With the density dependent bag constant, we have found that the maximum 
gravitational mass of SQS is about 1.62 M , while with the fixed bag constant, it is about 
1.33 M . 

We have plotted the gravitational mass of SQS as a function of the radius (mass-radius 
relation) for the magnetic fields B = and 5 x 10 18 G in Fig. [9j It is seen that for all cases 
of SQS, the gravitational mass increases by increasing the radius. In Fig. |9l we have also 
compared our results for the density dependent case of bag constant with those of density 
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independent case. We can see that for the case of fixed bag constant, the increasing rate 
of gravitational mass versus radius is higher than that of density dependent case. However, 
it will be more constructive to consider the effects of rotation on the properties of the star 
which is beyond our present investigation. Some authors have shown that considering the 
rot ation of the star leads to the larger maximum mass for strange quark stars (Shen et 



al. 



20051 ). 



IV. SUMMARY AND CONCLUSIONS 



We have investigated a cold static strange quark star in the presence of a strong magnetic 
field. For this purpose, some of the bulk properties of the strange quark matter such as 
the energy density and equation of state have been computed using the MIT bag model 
with the density dependent bag constant. Calculations of the energy for different magnetic 
polarization in the presence of a magnetic field demonstrated that as the density increases, 
the minimum point of energy shifts to the lower values of the polarization. We have shown 
that the value of the polarization parameter decreases by increasing the density, and it also 
increases by increasing the magnetic field. Our results at B = 5 x 10 18 G show that for 
both cases of the density dependent bag constant and fixed bag constant, the total energy 
density have an increasing rate by increasing the density. We have shown that there is 
a ferromagnetic phase transition at high magnetic fields. Our computations indicate that 
the pressure increases by increasing the density and magnetic field. In this work, we have 
also studied the structure properties of the strange quark stars. Our results show that 
the gravitational mass of the strange quark star increases by increasing the central energy 
density. It was shown that this gravitational mass reaches a limiting value at higher values 
of the central energy density. We have shown that the maximum mass of the strange quark 
star reduces by presenting the magnetic field. Finally, a comparison has been also made 
between the results of density dependent bag constant and those of fixed bag constant. Our 
calculation with the density dependent bag constant shows a higher maximum mass with 
respect to that of fixed bag constant. 

One of the possible astrophysical implications of our results is calculation of the surface 
redshift (z s ) of SQS. This parameter is of special interest in astrophysics and can be obtained 
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from the mass and radius of the star using the following relation (Camenzind 20071 ] ) 



2GM S i , s 

Z - = (1 - S^' 1 - L < 24) 

Our results corresponding to the maximum mass and radius of SQS calculated by density 
dependent bag constant lead to z s = 0.534 ms _1 for the magnetic field B = 5 x 10 18 G. 
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FIG. 1: Total energy density (etot) as a function of the polarization parameter (£) for different 
densities (p). 
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FIG. 3: Polarization parameter (£) corresponding to the minimum points of energy density versus 
the magnetic field (B) for different values of density (/?). 
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FIG. 4: Total energy density versus density (p) calculated by density dependent bag constant (full 
curve) at B = 5 x 10 18 G. The results for B = 90 ^ (dashed curve) at B = 5 x 10 18 G have 
been also given for comparison. 
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FIG. 5: Phase diagram for the strange quark matter in the presence of a strong magnetic field. 
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FIG. 7: The equation of state of SQM in the case of density dependent bag constant (full curve) 
at B = 5 x 10 18 G. The results for the case of fixed bag constant (B = 90 j^f) (dashed curve) at 
B = 5 x 10 18 G have been also given for comparison. 
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FIG. 8: Gravitational mass versus the central energy density (e c ) at B = (full curve) and 
B = 5x 10 18 G (dashed curve). The results for B = 90 ^ (dashed dotted curve) at B = 5x 10 18 G 
have been also given for comparison. 
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FIG. 9: The gravitational mass versus radius at B = (full curve) and B = 5 x 10 18 G (dashed 
curve). The results for B = 90 (dashed dotted curve) at B = 5 x 10 18 G have been also given 
for comparison. 
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